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Abstract 

Rich clusters of galaxies, the largest virialized systems known, 
place some of the most powerful constraints on cosmology. I discuss 
below the use of clusters of galaxies in addressing two fundamental 
questions: What is the mass-density of the universe? and how is the 
mass distributed? I show that several independent methods utilizing 
clusters of galaxies — cluster dynamics and mass-to-light ratio, baryon 
fractions in clusters, and cluster evolution — all indicate the same ro- 
bust result: the mass-density of the universe is low, f2 m ~0.2, and the 
mass approximately traces light on large scales. 



1 Introduction 

Theoretical arguments based on standard models of inflation, as well as on 
the demand of no "fine tuning" of cosmological parameters, predict a flat 
universe with the critical density needed to just halt its expansion. The crit- 
ical density, 1.9 x 10 _29 h 2 g cm -3 (where h refers to Hubble's constant, see 
below), is equivalent to ~10 protons per cubic meter; this density provides 
the gravitational pull needed to slow down the universal expansion and will 
eventually bring it to a halt. So far, however, only a small fraction of the 
critical density has been detected, even when all the unseen dark matter in 
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galaxy halos and clusters of galaxies is included. There is no reliable indi- 
cation so far that most of the matter needed to close the universe does in 
fact exist. Here we show that several independent observations of clusters 
of galaxies all indicate that the mass density of the universe is sub-critical. 
These observations include the mass and mass-to-light ratio of clusters and 
superclusters of galaxies, the high baryon fraction observed in clusters, and 
the evolution of the number density of massive clusters with time; the lat- 
ter method provides a powerful measure not only of the mass-density of the 
universe but also the amplitude of the mass fluctuations. The three indepen- 
dent methods- all simple and robust- yield consistent results of a low-density 
universe with mass approximately tracing light on large scales. 

2 Cluster Dynamics and the Mass-to-Light 
Ratio 

Rich clusters of galaxies are the most massive virialized objects known. Clus- 
ter masses can be directly and reliably determined using three independent 
methods: (i) the motion (velocity dispersion) of galaxies within clusters re- 
flect the dynamical cluster mass, within a given radius, assuming the clusters 
are in hydrostatic equilibrium [1-3]; (ii) the temperature of the hot intraclus- 
ter gas, like the galaxy motion, traces the cluster mass [4-6]; and (iii) grav- 
itational lensing distortions of background galaxies can be used to directly 
measure the intervening cluster mass that causes the distortions [7-10]. All 
three independent methods yield consistent cluster masses (typically within 
radii of ~1 Mpc), indicating that we can reliably determine cluster masses 
within the observed scatter (~ ± 30%). 

The simplest argument for a low density universe is based on summing 
up all the observed mass(associated with light to the largest possible scales) 
by utilizing the well-determined masses of clusters. The masses of rich clus- 
ters of galaxies range from ~10 14 to 10 15 h" 1 M within l^li" 1 Mpc radius 
of the cluster center (where h = H /100 km s _1 Mpc -1 denotes Hubble's 
constant). When normalized by the cluster luminosity, a median mass-to- 
light ratio of M/L B ~ 300 ± lOOh in solar units (M /L ) is observed for rich 
clusters, independent of the cluster luminosity, velocity dispersion, or other 
parameters [3,11]. (Lb is the total luminosity of the cluster in the blue band, 
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corrected for internal and Galactic absorption.) When integrated over the 
entire observed luminosity density of the universe, this mass-to-light ratio 
yields a mass density of p m ~ 0.4 x lCT 29 h 2 g cm" 3 , or a mass density ratio 
of Q m = p m / Peru — 0.2 ±0.1 (where p crit is the critical density needed to 
close the universe). The inferred density assumes that all galaxies exhibit 
the same high M/Lb ratio as clusters, and that mass follows light on large 
scales. Thus, even if all galaxies have as much mass per unit luminosity 
as do massive clusters, the total mass of the universe is only ~20% of the 
critical density. If one insists on esthetic grounds that the universe has a 
critical density (Q m = 1), then most of the mass of the universe has to be 
unassociated with galaxies (i.e., with light). On large scales (>1.5 h" 1 Mpc) 
the mass has to reside in "voids" where there is no light. This would imply, 
for Q m = 1, a large bias in the distribution of mass versus light, with mass 
distributed considerably more diffusely than light. 

Is there a strong bias in the universe, with most of the dark matter 
residing on large scales, well beyond galaxies and clusters? An analysis of 
the mass-to-light ratio of galaxies, groups, and clusters by Bahcall, Lubin 
and Dorman [11] suggests that there is not a large bias. The study shows 
that the M/Lb ratio of galaxies increases with scale up to radii of R ~0.2 
h" 1 Mpc, due to very large dark halos around galaxies [see also 12,13]. The 
M/L ratio, however, appears to flatten and remain approximately constant 
for groups and rich clusters from scales of ~0.2 to at least 1.5h _1 Mpc and 
possibly even beyond (Figure 1). The flattening occurs at M/L B ^ 200-300h, 
corresponding to fl m ~ 0.2. (An M/L B c± 1350 h is needed for a critical 
density universe, Q m =l.) This observation contradicts the classical belief 
that the relative amount of dark matter increases continuously with scale, 
possibly reaching Q m =l on large scales. The available data suggest that most 
of the dark matter may be associated with very large dark halos of galaxies 
and that clusters do not contain a substantial amount of additional dark 
matter, other than that associated with (or torn-off from) the galaxy halos, 
plus the hot intracluster gas. This flattening of M/L with scale, if confirmed 
by further larger-scale observations, suggests that the relative amount of dark 
matter does not increase significantly with scale above ~0.2h _1 Mpc. In that 
case, the mass density of the universe is low, fi m ~0.2-0.3, with no significant 
bias (i.e., mass approximately following light on large scales). 

Recently the mass and mass-to-light ratio of a supercluster of galaxies, 
on a scale of ~6h _1 Mpc, was directly measured using observations of weak 
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gravitational lensing distortion of background galaxies (Kaiser et al. [14]). 
The results yield a supercluster mass-to-light ratio (on 6h _1 Mpc scale) of 
M/L B = 280 ±40h, comparable to the mean value observed for the three indi- 
vidual clusters that are members of this supercluster. These results provide 
a powerful confirmation of the suggested flattening of M/Lb (R) presented 
in Figure 1 (Bahcall et al. [11,15]). The recent results confirm that no sig- 
nificant amount of additional dark matter exists on large scales. The results 
also provide a clear illustration that mass approximately traces light on large 
scales and that Q m is low, as suggested by Figure 1. 

3 Baryons in Clusters 

Clusters contain many baryons, observed as gas and stars. Within 1.5h -1 
Mpc of a rich cluster, the X-ray emitting gas contributes ~6h _L5 % of the 
cluster virial mass [16-18]. Stars contribute another ~2%. The baryon frac- 
tion observed in clusters is thus: 

tt h /n m > 0.06h- 15 + 0.02 (1) 

The observed value represents a lower-limit to the baryon fraction since we 
count only the known baryons observed in gas and stars; additional baryonic 
matter may of course exist in the clusters. Standard Big Bang nucleosynthesis 
limits the baryon density of the universe to [19,20]: 

Vl b ~ 0.015h" 2 (2) 

These facts suggest that the baryon fraction observed in rich clusters (eq. 1) 
exceeds that of an f2 m =l universe (f2 6 /(f2 m =l)~0.015h~ 2 ; eq. 2) by a factor 
of >3 (for h>0.5). Since detailed hydrodynamic simulations [16,18] show 
that baryons do not segregate into rich clusters, the above results imply that 
either the mean density of the universe is lower than the critical density by 
a factor of >3, or that the baryon density is much larger than predicted by 
nucleosynthesis. The observed high baryonic mass fraction in clusters (eq. 1), 
combined with the nucleosynthesis limit (eq. 2), suggest (for h~0.65±0.1): 

On < 0.3. (3) 

This upper limit on Q m is a simple, model-independent and thus powerful 
constraint: a critical density universe is inconsistent with the high baryon 
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fraction observed in clusters (assuming the nucleosynthesis limit on f^, eq. 
2). [A universe dominated by hot dark matter that is too hot to clump into 
clusters may survive this critical test; however, such a hot universe is unable 
to form the observed galaxies at high redshift, and is therefore unlikely] 

4 Evolution of Cluster Abundance 

The observed present-day abundance of rich clusters of galaxies places a 
strong constraint on cosmology: o" 8 f2 m °- 5 ~0.5, where <7 8 is the rms mass 
fluctuations on 8h _1 Mpc scale, and Q, m is the present cosmological density 
parameter [21-26]. This constraint is degenerate in Q m and a s ] models with 
Q m =l, o" 8 ~0.5 are indistinguishable from models with f2 m ~0.25, <x 8 ~l. (A 
og^l universe is unbiased, with mass following light on large scales since 
galaxies (light) exhibits cr 8 (galaxies) ~1; <r 8 ~0.5 implies a biased universe 
with mass distributed more diffusely than light). 

The evolution of cluster abundance with redshift, especially for massive 
clusters, breaks the degeneracy between Q m and a 8 [see, e.g. 23,24,27-34]. 
The evolution of high mass clusters is strong in Q m =l, low-<r 8 (biased) Gaus- 
sian models, where only a very low cluster abundance is expected at z > 0.5. 
Conversely, the evolution rate in low-f2 m high-a 8 models is mild and the 
cluster abundance at z > 0.5 is much higher than in Q m =l models. 

In low-density models, density fluctuations evolve and freeze out at early 
times, thus producing only relatively little evolution at recent times {z < 1). 
In an Vt m =l universe, the fluctuations start growing more recently thereby 
producing strong evolution in recent times; a large increase in the abundance 
of massive clusters is expected from z ~ 1 to 2; ~ 0. In a recent study by 
Bahcall, Fan, and Cen [31] we show that the evolution is so strong in f2 m =l 
models that finding even a few Coma-like clusters at z > 0.5 over ~10 3 deg 2 
of sky contradicts an Q m =l model where only ~10~ 2 such clusters would 
be expected (when normalized to the present-day cluster abundance). The 
evolution of the number density of Coma-like clusters was determined from 
observations using the CNOC cluster sample to z < 0.5 and compared with 
cosmological simulations [30-32]. The data show only a slow evolution of 
the cluster abundance to z ~ 0.5, with ~10 2 times more clusters observed at 
these redshifts than expected for Q m —1. The results yield f2 m ~0.3±0.1. 

The evolutionary effects increase with cluster mass and with redshift. The 
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existence of the three most massive clusters observed so far at z ~ 0.5 — 0.9 
places the strongest constraint yet on Q m and as- These clusters (MS0016+16 
at z = 0.55, MS0451-03 at z = 0.54, and MS1054-03 at z = 0.83, from the 
Extended Medium Sensitivity Survey, EMSS [35,36]), are nearly twice as 
massive as the Coma cluster, and have reliably measured masses (including 
gravitational lensing masses, temperatures, and velocity dispersions; [34,37- 
40]. These clusters posses the highest masses (>8 xlO 14 h -1 M within 1.5 
h _1 comoving Mpc radius), the highest velocity dispersions (>1200 km s _1 ), 
and the highest temperatures (>8 keV) in the z > 0.5 EMSS survey. The 
existence of these three massive distant clusters, even just the existence of 
the single observed cluster at z — 0.83, rules out Gaussian Vt m =l models 
for which only ~10~ 5 z ~ 0.8 clusters are expected instead of the 1 cluster 
observed (or ~10~ 3 z > 0.5 clusters expected instead of the 3 observed). (See 
Bahcall & Fan [34]). 

In Figure 2 we compare the observed versus expected evolution of the 
number density of such massive clusters. The expected evolution is based 
on the Press-Schechter [41] formalism; it is presented for different fl m val- 
ues (each with the appropriate normalization og that satisfies the observed 
present-day cluster abundance, <x 8 ~0.5 f2 m ~ a5 ; [23,26]). The model curves 
range from fl m =0.1 (crg~1.7) at the top of the figure (flattest, nearly no 
evolution) to Q m =l ((X8~0.5) at the bottom (steepest, strongest evolution). 
The difference between high and low Q m models is dramatic for these high 
mass clusters: Vl m =l models predict ~10 5 times less clusters at z ~ 0.8 than 
do f2 m ~0.2 models. The large magnitude of the effect is due to the fact 
that these are very massive clusters, on the exponential tail of the cluster 
mass function; they are rare events and the evolution of their number den- 
sity depends exponentially on their "rarity", i.e., depends exponentially on 
a 8 2 oc VL m [32,34]. The number of clusters observed at z ~ 0.8 is consistent 
with fi m ~0.2, and is highly inconsistent with the ~10~ 5 clusters expected 
if Q m =l. The data exhibit only a slow, relatively flat evolution; this is ex- 
pected only in low-f2 m models. Vt m =l models have a ~10~ 5 probability of 
producing the one observed cluster at z ~ 0.8, and, independently, a ~10~ 6 
probability of producing the two observed clusters at z ~ 0.55. These results 
rule out f2 m =l Gaussian models at a very high confidence level. The results 
are similar for models with or without a cosmological constant. The data 
provide powerful constraints on Q m and cr 8 : Q m =0.2~^Q'\ 5 and cr 8 = 1.2 ± 0.3 
(68% confidence level) [34]. The high <r 8 value for the mean mass fluctuations 
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indicates a nearly unbiased universe, with mass approximately tracing light 
on large scales. This conclusion is consistent with the suggested flattening of 
the observed M/L ratio on large scales (Figure 1, section 2). 

The rate of evolution of the cluster abundance depends strongly on <7g, 
as was shown by Fan et al. [32]: dlogn/dz oc-a^ 2 ; low-cr 8 models evolve 
exponentially faster than high-cr 8 models, for a given mass cluster. The 
strong exponential dependence arises because clusters of a given mass repre- 
sent rarer density peaks in low-as models. We show [32] that the evolution 
rate at z<l is relatively insensitive to the density parameter Q m or to the 
exact shape of the power spectrum. This is illustrated in Figure 3, where 
we present the cluster abundance ratio, n(z~0.8)/n(z~0), as a function of 
0"s for all Q m values. The strong exponential evolution rate for low-<r 8 is 
clearly distinguished from the nearly no evolution expected for a 8 ~l (f° r 
any Q m ). The dependence on Q m is only secondary. This method thus pro- 
vides a unique tool for determining <r 8 — by simply measuring the observed 
ratio of cluster abundance at different redshifts. When combined with the 
observed normalization of the present-day cluster abundance, <7 8 f2 m a5 ~ 0.5, 
the density parameter Q m can be determined. The observed evolution rate 
implied by the existence of the massive high redshift clusters discussed above 
is shown in Figure 3 — it indicates only minimal evolution thus a high <r 8 
value of ~1. 

In Figure 4 we summarize the above results for Q m and <r 8 as determined 
from cluster evolution and compare the results with the other two indepen- 
dent methods discussed in sections 2 and 3. We show in Figure 4 the Vt m — a 8 
band allowed by the present-day cluster abundance f2 m °' 5 ~0.5 erg -1 , and the 
band allowed by the existence of the high- redshift massive clusters [34] . The 
intersection of these two bands provides the cluster evolution constraints dis- 
cussed above; the only allowed range for Gaussian models is low-f2 m , <t 8 ~1. 
We also present the Vt m constraints determined from the two entirely inde- 
pendent methods discussed in §§2-3: the high baryon fraction observed in 
clusters, which yields f2 m <0.3 (assuming the nucleosynthesis limit for the 
baryon density); and cluster dynamics which yields Vt m ~0.2. All these in- 
dependent methods yield consistent results: fi m ~0.2 ±0.1 and (j 8 ~1.2±0.2 
(1 a level). f2 m =l models are highly incompatible with these results. 
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5 Summary 



We have shown that several independent observations of clusters of galaxies 
all indicate that the mass-density of the universe is sub-critical: fi m ~0.2±0.1. 
The results also suggest that mass approximately traces light on large scales. 
A summary of the results, presented in Figure 4, is highlighted below. 

1. The mass-to-light ratio of clusters (and superclusters) of galaxies and 
the suggested flattening of the mass-to-light ratio on large scales suggest 
Q m ~0.2±0.1. 

2. The high baryon fraction observed in clusters of galaxies suggests 

n m <o.3. 

3. The weak evolution of the observed cluster abundance to z ~ 1 pro- 
vides a robust estimate of f2 m ~0.2^Q;} 5 , valid for any Gaussian models. An 
Q m =l Gaussian universe is ruled out as a <1CT 6 probability by the cluster 
evolution results (Figure 2-4). 

4. All the above-described independent measures are consistent with 
each other and indicate a low-density universe with fi m ~0.2±0.1 (Figure 
4). f2 m =l models are ruled out by the data. While non-Gaussian initial 
fluctuations, if they exist, will affect the cluster evolution results, they will 
not affect arguments (1) and (2) above. Gaussian low-density models (with or 
without a cosmological constant) can consistently explain all the independent 
observations presented here. These independent cluster observations indicate 
that we live in a lightweight universe with only ~20%-30% of the critical 
density. Thus, the universe may expand forever. 
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Figure 1: The dependence of mass-to-light ratio, M/Lb, on scale, R, for average spiral 
galaxies (spiral symbols), elliptical galaxies (elliptical symbols), and groups and clusters 
(filled circles). Adapted from Bahcall, Lubin and Dorman [11]; [15]. The large scale 
point at ~15 h _1 Mpc represents Virgo cluster infall motion results [11]. The location 
of TO = 1 and f2 TO = 0.3 are indicated by the horizontal lines. A flattening of M/Lb is 
suggested at il m ~ 0.2 ± 0.1. A recent result for a supercluster using weak gravitational 
lensing (R~6h _1 Mpc; [14]), M/LB=280±40h, is consistent with the suggested flattening 
of M/L(R). 
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Figure 2: Evolution of the number density of massive clusters as a function of redshift: 
observed versus expected (for clusters with mass >8 x 10 14 h _1 Mq within a comoving 
radius of 1.5h _1 Mpc). (Adapted from Bahcall and Fan 1998 [34]. The expected evolu- 
tion is presented for different TO values by the different curves. The observational data 
points (see text) show only a slow evolution in the cluster abundance, consistent with tt m 
c±0.2~^q \ 5 . Models with O m =l predict ~10 5 fewer clusters than observed at z ~0.8, and 
~10 3 fewer clusters than observed at z ~0.6. 
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n=o. 2-1.0 (n + a = 1 ) 




Figure 3: Cluster abundance ratio, n(z — 0.8)/n(z = 0), vs. as from Press-Schcchtcr 
(solid curve, for mean of all Sl m 's) for clusters with mass >8 x 10 14 h _1 M . Filled circles 
represent f2 m 's from 0.2 to 1 (bottom to top). (Dashed curves represent the mass threshold 
range of 7 to 10 x 10 14 M Q , top and bottom, respectively). The data (Figure 2, §4) are 
shown by the shaded region (68% level). Similar results are obtained for A = 0. 
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Figure 4: Constraining the mass-density parameter, Q m , and the mass fluctuations on 
8h _1 Mpc scale, as, from several independent observations of clusters: cluster dynamics; 
baryon fraction in clusters; present-day cluster abundance (z ~ 0); and cluster abundance 
at redshift z ~ 0.7. (The latter two abundances yield the cluster evolution constraints 
shown in Figure 2; see text). All these model-independent observations converge at the 
allowed range of Q m = 0.2 ± 0.1 and cr 8 = 1.2 ± 0.2 ( 68% confidence level). The dotted 
lines illustrate the mean microwave fluctuations constraints, based on the COBE satellite 
results, for a Cold-Dark-Matter model with h = 0.7 (with and without a cosmological 
constant, denoted as LCDM and OCDM respectively. Both models are consistent, within 
their uncertainties, with the best-fit fl m — as regime of the cluster observations). 
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